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Abstract 

The symmetry group of the staggered Fermion transfer matrix in 
a spatial direction is constructed at finite temperature. Hadron-like 
operators carrying irreducible representations of this group are written 
down from the breaking of the zero temperature group. Analysis of 
the correlators in a free fermion theory suggests new measurements 
which can test current interpretations. 
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1 Introduction 



Lattice simulations of field theories in equilibrium at finite temperature (T) 
use a discretisation of the Euclidean formulation for partition functions — 



Z((3) = J Ztyexp - J dtj d 3 x£(<p) 



(1) 



where is a generic field, L the Lagrangian density, and the Euclidean "time" 
runs from to (3 = 1/T. The path integral is over Bosonic (Fermionic) field 
configurations which are periodic (anti-periodic) in Euclidean time. Due to 
the lack of symmetry between the space and Euclidean time directions in 
eq. (|I|), this problem has only a subgroup of the full 4-dimensional rotational 
symmetry of the T = Euclidean theory. In this paper we focus on the lattice 
discretised problem, where all continuum symmetries break to a discrete 
subgroup. 

It is possible to write the partition function of eq. ([!]) as the trace of 
the transfer matrix in one of the spatial directions. The symmetry groups 
we examine leave such a transfer matrix invariant. The eigenvectors of the 
matrix carry irreducible representations (irreps) of these symmetry groups. 

Thermodynamics depends only on the leading eigenvalue, which always 
belongs to a scalar representation of the symmetry group. Hence the group 
theory is not crucial for the study of properties such as the phase structure, 
transition temperature, T c , and other thermodynamic quantities. In fact, 
extensive measurements have been made of T c for pure gauge theories, and 
those with massless Fermions [|IJ], and our group theoretical analysis adds 
very little to this. 

However, the symmetry properties are crucial to the study of screening 
correlation functions and the determination of screening masses. These can 
be written in terms of the ratio of the largest and an appropriate other eigen- 
value of the transfer matrix. The significance of the equality (or otherwise) of 
two screening masses will depend on whether or not the correlation function 
lies in the same irrep of the symmetry group of the transfer matrix. In the 
gauge sector of the theory this analysis has been carried out and applied to 
the study of screening masses [|2j , to demonstrate dimensional reduction in a 
fully non-perturbative manner. 

Screening masses obtained from correlation functions built out of stag- 
gered Fermion field operators have also been extensively studied in the past 
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[||, [|, [|. Screening masses in the high temperature phase of QCD seem to 
approach those expected from free field theory as T — > 2T C @, [|, |[]. Some 
other measurements which seem to indicate that the picture may be more 
complicated ||, also turn out to be explained in terms of weakly interacting 
quarks f7j. All these studies have relied entirely on the T = analysis of the 
lattice symmetry group of staggered Fermions. 

In this paper we present the first analysis of the symmetries of the corre- 
sponding finite temperature problem. We find that all the screening masses 
measured till now see only one of the representations of the symmetry group. 
Many other masses can be studied, and are likely to yield further informa- 
tion about the theory. The free field theory of these other representations is 
worked out. 

One observation arising from the application of these group theoretical 
results to previous simulations is worth mentioning in the introduction. Since 
the T = scalar and pseudo-scalar mesons, and the symmetric linear com- 
binations of the three components of the vector and pseudo- vector lie in the 
same irreducible representation (irrep) of the point group of a T > spatial 
slice through the lattice, they must have degenerate masses in a free fermion 
theory. When interactions switch on, the relevant symmetry becomes that 
of an enveloping group, and the four degenerate masses split into two pairs 
of degenerate masses. Observation of such a splitting for T < 2T C |3|, [|, [|. || 
must then be interpreted as evidence for interactions 0. Nothing further can 
be said purely from the study of these correlators. Whether the spectrum of 
screening masses comes from a weakly interacting effective theory, or whether 
it is very similar to the spectrum at zero temperature, are questions which 
can only be answered by measuring the masses in the other representations 
which we write down explicitly. 

In Section we present a brief review of the symmetries of staggered 
Fermions at T = 0. This serves to set up the notation, and indicates what 
changes to expect at finite temperature. Section contains our main results 
on the characterisation of the group of symmetries of the spatial transfer 
matrix at T > and its irreducible representations (irreps). Free field theory 
results for the screening masses is discussed in Section f|, where presently 
available data is also discussed. Two appendices contain the technical details 
of induced representations and character tables for the irreps of mesons. 
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2 Symmetries of Staggered Fermions 



In this section we review the breaking of continuum spin-flavour symmetries 
for lattice staggered Fermions f| at zero temperature, and identify how this 
pattern changes at finite temperature. The continuum symmetry for four 
flavours of Fermions is SU r (2) <g> Uf(4), where the first factor is the rotational 
symmetry, and the second is the flavour symmetry. We follow the notational 
conventions of || 10 



At T = we are interested in the symmetries of Fermion operators which 
have zero momentum in the directions orthogonal to the Euclidean time t — 



Xa 



y x 



nil 



X (r)T^T™>T™* = $>(x + 2ma). (2) 



Here the index A denotes the corners of the hypercube on which the ap- 
propriate component of the quark field resides, a is the lattice spacing, Tj 
are the generators of translations in the i-th direction, and we have assumed 
that there are periodic boundary conditions in all directions on the slice. In 
writing eq. (0), we have chosen to study correlation functions of operators 
separated in the time direction. Due to the 4-dimensional discrete rotational 
symmetry of the T = theory, we could have chosen to study propagation 
in any other direction with the same result. 



Operation 



Action 



FT 



R 



S, 



c 



K(x) 



x 



= \ 1 ± ri K (x)ri x (x) =F + 

+r} K (x)r]x(x)( K (x)(x(x) 

X(x) -> rj 4 (x)x(Ix) 
X(x) -> C^(x)x(x + dp) 
X(x) -> e(x)x(x) 



X(x) 



X{x) 



Table 1: Symmetry operations on staggered Fermions. The upper (lower) 
signs in TZ are used when k > A (k < A). Here e(x) = (— i) x i+ a: 2+a;3+^4 ) 
rii{x) = n fe<i (-l)^ and Q(x) = n^-l)**- 
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The symmetry elements of the theory are listed in Table |I[ For stag- 
gered Fermions, the shifts by one lattice spacing, S^, are mixed flavour and 
translation operations. Pure translations are = S 2 . We have chosen the 
transfer matrix T to be T t . Nothing would have changed, at T = 0, if we 
had instead chosen T to be T z . 

Discrete flavour operations, = S^T^ 1 / 2 , are vectors under rotations, 
generated by R K \ and transform as 

Ri/^kRij = 8ik^j + 8jk^i + \tijk\^k- (3) 

Here and elsewhere, Greek indices run from 1 to 4; Latin indices over the 
three directions summed in eq. (0) or its analogue. A subgroup Ub{1) of 
the continuum flavour group remains unbroken on the lattice; this charge 
corresponds to the Fermion number q. The representations of S M in an irrep 
with Fermion number q, T> g (H M ), obey the relation 

V q (Z,)V q (~ u ) = ^V q (E u )V q (E,). (4) 

Inversion, /, commutes with S 4 , and anti-commutes (commutes) with the 
other 2^ in representations with odd (even) values of q. Parity is defined by 
P = S 4 J. The remaining discrete symmetry is that of charge-conjugation, 
C. 

The symmetries of T are the rest-frame group 

RF(E^R kh I,C)®U B (l). (5) 

We have used the notation G(X) to mean the group G generated by the oper- 
ation^) X. A subgroup of RF is the group of isometries of the lattice, called 
the geometric rest frame group, GRF (S M , Rm, I). In turn, GRF contains 
the time slice group, which is the point group of the lattice — 

TS (R kl , I) = O h (R kl , I) = (R kl ) ® Z 2 (I). (6) 

This chain of groups builds up to the continuum symmetry group — 

TS C GRF c RF c SU d {2) <g> U B (1), (7) 

where SUd{2) is the diagonal subgroup of the direct product SU r {2)®SUf(2) 
of rotations and flavour. The breaking of SUf(4) to S77/(2) is specified 
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by requiring that the fundamental of SUf(4) break into the irrep (|, \) of 
5*7,(2). 

All correlation functions block diagonalise into irreps of TS. This group, 
Oh, is the group of symmetries of a cube. It has 48 elements in 10 conju- 
gacy classes [13[]. It has four one-dimensional irreps Af and Af, two two- 



dimensional irreps and four three-dimensional irreps F x and F 2 . The 
physical interpretation of each mass is obtained by tracing the descent of 
the irrep of TS through the whole chain in eq. (|7|) from the irreps of the 
continuum symmetry 5*0(4) g) SU(4:). This is done in |9|, p70f| . See also |T 
for some details of the treatment of correlation functions. 

For the study of equilibrium finite temperature, T > 0, physics we are 
interested in screening masses and screening correlation functions, i.e., in the 
eigenvalues of the transfer matrix in spatial directions. Two distinctions from 
the T = theory should be borne in mind. 

The first is that there are anti-periodic boundary conditions in the Eu- 
clidean time direction on Fermions. As a result the lowest Fourier component 
has a non-vanishing momentum in this direction — 

Tf Nt/2 X (x)T t Nt/2 = e*" X {x), T t ~ Nt/2 X (x)T t Nt/2 = e^ X (x), (8) 

where N t is the number of lattice points in the time direction. This is a trivial 
change. For Fermion bilinear operators it makes no difference. Operators 
with an odd number of Fermion fields are treated slightly differently. For 
example, the projection on the lowest momentum state of a Fermion field is 
not written as in eq. (Q), but as 

XA = Ee 2 ™ tM x(x + 2ma), (9) 

m 

where m runs over the coordinates in a spatial slice, i.e., over two spatial 
directions and the temporal direction [|]. The phase factor in the sum is just 
the statement that the lowest Matsubara frequency for Fermions is ttT. 

In this paper we shall concern ourselves with the second, and more im- 
portant, difference — the isometries of a slice of the lattice. Since we are 
interested in screening masses, we consider slices through the lattice orthog- 
onal to one of the spatial directions, say the ^-direction (as in eq. || above). 
Then the isometries of the z-slice generate 

TS = D h 4 = D 4 (R xy , R 2 xt ) ® Z 2 (I). (10) 
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The identification of this group is easy, because it differs from Oh (eq. ||) by 
the fact that rotations of it/2 in the xt and yt planes is not allowed. The 
spectrum of the screening masses requires a classification by the irreps of D\. 
The continuum symmetry will be built up by the group chain 

TS c GRF c RF c C ® U B {1) c SU d {2) ® U B {1), (11) 

where C = 0(2) <g> Z 2 (I) is the invariance group of a cylinder. To generate 
each of the lattice groups in the chain, we use the construction at T = 0, 
only leaving out odd powers of the rotations R^ t . 

D\ has 16 elements in six conjugacy classes ||13|| . There are eight one- 



dimensional irreps labelled Af, A 2 , Bf and B 2 -, and two two-dimensional 
irreps E^. The reductions of the irreps of Oh to D\ is as — 

Af — > A p , A P — > B[, 

F[ ^ A!-®E P , F 2 P - B* @E P , 

E p -> Af ©Sf. (12) 

More details can be found in 0. 

In the rest of this paper we shall give these decompositions of meson 
and hadron operators using the language of the T = theory. This calls 
for some care in the interpretation of results — although we shall talk of 
charge conjugation, C, and parity, P, and the operators will have the same 
structure and algebra as in the T = theory, they may represent quite 



different physical quantities [12 



3 The Symmetry Group at T > 0. 

In this section the symmetry groups are written down. The representation 
theory of these groups in the meson (quark-antiquark) sector is examined in 
detail. The symmetries of the quark fields are also examined briefly, and the 
representation theory in the baryon sector is dealt with in less detail. 

3.1 Meson Operators 

In a meson representation, the quark number q = 0. As a result, the repre- 
sentants of the flavour generators commute. Consequently C and / 
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G 





G 


D 4 


Meson 


1 


A x 


lo 


Ax 


Yes 


V 


A 2 


li 


Bi 




2 


E 


lo + li 


A! + B x 




3 


Fi 


2 2 + 1 6 


A 2 + E 


Yes 


3' 


F 2 


2 2 + 1 7 


B 2 + E 




3" 


Fx 


2 + l 2 


A 2 + E 


Yes 


3"' 


F 2 


2 + l 3 


B 2 + E 




3"" 


Ax + E 


2 4 + 1 4 


2A 1 + B x 


Yes 


3""' 


A 2 + E 


2 4 + l 5 


A 1 + 2Bt 




6 


Fx + F 2 


2 a + 2 3 + 2 5 


A 2 + B 2 + 2E 


Yes 



Table 2: Irreps of G, denned in eq. (p~4|) , and their reduction at finite tem- 
perature to irreps of G, defined in eq. flT5|). The irreps which are realised for 
mesons are marked. Meson states do not exhaust all the irreps of G, O or 
G, but do exhaust all the irreps of D4. 



commute with X k . At T = 0, it has been shown that [ID 

RF = GRF(X„,R kl ,I)<g,Z 2 (C), (13) 

where 

GRF (X^ R kh I) = G (X kl R kl ) ® Z 2 (I) ® Z 2 (X 1 X 2 X 3 ) ® Z 2 (X 4 ), (14) 

with X k = X fc XiX 2 X 3 . The irreps of GRF are denoted r a4ai23 , where r 
denotes an irrep of G, and <r 4 and <7i 2 3 are signs which denote the irreps of 
the Z 2 factor groups generated by X 4 and XiX 2 X 3 respectively. 

Next we identify the group G. The X k generate a 4 element Abelian 
group called the Viergruppe, V = Z 2 ® Z 2 . This is a normal subgroup 
of G. The transformation properties of X k under rotations, eq. (|3p, show 
that G is the semi-direct product G = V(X k ) tx O. Now, the cubic group 
O = V(Rl l ) M 5*3, where the normal subgroup, V(R kl ) is generated by the 
three rotations by angle tc ||14|| , and the other factor is the permutation group 
of 3 elements. From eq. (|3|) it is clear that V{R\{) has trivial action on V(X k ), 
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and we can write 



G(X k ,R kl ) = ( V{X k ) ® V{R kl ) I cxi S; 



(15) 



Since the normal subgroup is Abelian, the irreps of G can be efficiently 
generated by the method of induced representations. Details are given in 
Appendix KL where we recover the results of ITU . 



GRF 


G 




Operator 


1++ 

1+- 
3 //«+- 

3W/++ 


lo 
lo 
I4 

2 4 

I4 

2 4 


At 
At 
At 
At 
Bt 
At 
At 
Bt 


ExV4(x)(a(x)x(x)x(x) 
Ex ^(x)r] 3 (x)C 3 (x)x(x)x(x) 
Ex^(x)[vi(x)Ci(x) + ri2(x)( 2 (x)]x(x)x(x) 
Eaie(z)[»7i(a:)Ci(x) - V2(x)( 2 (x)}x(x)x(x) 
J2x^(x)r] 4 (x)C4(x)r] 3 (x)C3(x)x(x)x(x) 
T.x^{x)r] A (x)CA{x)[r]i{x)Ci{x) + rj 2 (x)C 2 (x)]x(x)x(x) 
Y^e(x)7i4,(x)C4(x)[Tii(x)Ci(x) - V2(x)( 2 (x)]x(x)x(x) 



Table 3: Representations of local staggered mesons. Only the At operators 
have been used in simulations till now. Reduction of three-link separated 
mesons follows an identical pattern and generates the opposite parity irreps 
of D\. 

This method makes it easy to construct the T > group, 

GRF = G®Z 2 {I)® Z 2 {X 1 X 2 X 3 ) ® Z 2 (X A ). (16) 

The rotation generators in G are R\ 2 and Rf 3 , and they generate the group 
D A . Since D A = V(Rl l ) ix Z 2 (R 12 ), we have 

G = (v(X k ) ® V(R 2 kl )j 1x1 Z 2 (R 12 ) (17) 

G has 32 elements in 14 conjugacy classes. The irreps can be constructed 
by the method of induced representations (see Appendix 0). There are 8 
one-dimensional and 6 two-dimensional irreps of G. 
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GRF 


G 




Operator 


3 + 


1 6 


^ 2 " 




Ex V3(x)x(x)D 3 x(x) 




2 2 


£T 




Ex ?7i,2(^)x(^)-Di, 2 x(^) 


3 


1 6 


A 2 




Ex ^(a;)C4(a;)^3(a;)x(^)^ , 3X(^) 




2 2 


E~ 




Ex V4(x)C4(x)rj 1:2 (x)x(x)D 1: 2x(x) 


3" 


1 2 


A 2 




Ex t(x)( 3 (x)x(x)D 3 x(x) 




2 


E- 




Ex e(^)Ci,2(^)x(^)-Di, 2 x(^) 


3 //-+ 


1 2 


A 2 




Ex Vi(x)C4(x)e(x)C 3 (x)x(x)D 3 x(x) 




2 


E~ 




Ex ?74(^)C4(^)e(^)Ci,2(a:)x(^)-Di, 2 x(^) 


6 


2 5 


A 2 




E x e(a:)[77i(x) + r} 2 (x)]rj 3 (x)x(x)D 3 x(x) 






B 2 




Ex e ( x )[Vi(x) - V2(x)]r] 3 (x)x(x)D 3 x(x) 




2i 


E~ 




J2 x e(x)r] 1 (x)rj 2 (x)x(x)D 2 x(x) : 1 <-> 2 




2 3 


E~ 




Exe(^)%(^)?7i(^)x(^)-Dix(^), 1 -> 2 


6 + 


2 5 


A 2 




(x)77 4 (a:)C4(a;)[77i(x) + ^(^^(^X^Ao:^) 






B 2 


Ex£ 


(x)77 4 (a:)C4(a:)[?7i(a:) - ^(^^(x^O^AjX^) 




2i 


E~ 


Ex 


e(x)r] 4 (x)C4(x)r] 1 (x)r] 2 (x)x(x)D 2 x(x), 1 <-> 2 




2 3 


E~ 


Ex 


e(x)r] 4 (a:)C4(a;)^3(^)^i(^)x(^)^ , iX(^) ) 1 -»■ 2 



Table 4: Representations of one-link separated staggered mesons. Here 
D^0(x) = 0(x + /2) + 0(x — /}). Reduction of two-link separated mesons 
follows an identical pattern and generates the opposite parity irreps of D\. 
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The content of the various GRF irreps is shown in Table 0. The reduction 
of irreps of G to those of O, G and D4 are performed using the character 
tables in Appendix |FJ. These reductions are consistent with those given in 
eq. (|T2D . The irreps obtained for mesonic (quark-antiquark) operators can 
be identified through the Clebsch-Gordan series for the GRF of Fermionic 
representations. 

We demonstrate the reduction of the l +± and 3"" +± irreps of GRF to 
the irreps of D\ in Table [| using the full set of local meson operators. All 
the operators which have been used to date for computing screening masses 
belong to the D\ irrep Af, and conversely, all the Af operators in Table |3] 
have been used in measurements. Notice, however, that this one irrep of TS 
descends from different irreps of GRF . The two Af irreps descending from 
1 must give degenerate masses^, as must the pairs descending from the 1 4 
and the 2 4 . However, there is no group theoretical necessity for the three 
pairs to have the same mass. 

The reduction of 3 _± , 3 //_± and 6 _± irreps obtained for one-link sepa- 
rated meson operators is given in Table |j. Reductions of two-link separated 
meson operators can also be read off from the structure of these reductions. 
The latter give positive I parity irreps of D\. Combining these two sets we 
have a set of Af, Bf and E ± irreps. Three link separated operators reduce 
in the same way as the local meson operators but give the opposite / parity. 
These two sets together give us the remaining irreps of D\, i.e., Af and Bf. 



3.2 Quark and Baryon Operators 

Quark and baryon operators carry odd Fermion charge. The representants 
of S M anti-commute and generate the 32 element Clifford group CL(A). Its 
commutator subgroup is isomorphic to Z 2 and its Abelisation, CL(4)/Z2 
is precisely the group V(X k ) <g> Z 2 (X 4 ) ® Z 2 (X 1 X 2 X 3 ) encountered as the 
normal subgroup of the mesonic rot at ion- shift group [jTJJ. Apart from the 



16 one-dimensional irreps of this group, CL(4) also has a four- dimensional 
quaternionic irrep familiar to us from the algebra of the Dirac matrices. 
The rot at ion- shift group for Fermionic operators at T = is 

G F = CL{±) M O = (CL(4) <g> V{R 2 kl )) ixa S 3 . (18) 



2 This is the phenomenon of "parity doubling" at high temperature. 
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GRF 


o h 




Component 


8 


At 


4" 


x(o) 




Al 


Al 


X(x + y + z) 






At 


X(x + v) 






E + 


{x(x + 1), -x(y + z)} 






A2 








E~ 





Table 5: Representations of the quark field xa- The "zero-momentum" pro- 
jection is performed as shown in eq. (|9|). The 8 of Gf is the same as the 80 
of the G F 

This group has 45 conjugacy classes, and hence 45 irreps. 40 of these have 
been identified as the irreps r°" 4<Jl23 of the mesonic rotation-shift group. The 
remaining 5 are obtained by inducing with the remaining non-trivial irrep of 
CL(4). This gives the five real irreps 8, 8', 16, 24 and 24'. 

The defining representation of Gf, 8, is given by zero momentum stag- 
gered Fermions on a time slice eq. (|2]). Under Oh the octet breaks || 

8 -> At + Ai+F+ + F^. (19) 

The Clebsch-Gordans for 8x8x8 show that only the 8, 8' and 16 are found 
as irreps of baryons |§. The local baryon operators built from staggered 
Fermions transform as the Af component of the 8. 

The rot at ion- shift group for Fermionic operators at T > is 

G F = CL(4) M D 4 = (CL(4) ® V(R 2 kl )) M Z 2 . (20) 

This group has 61 conjugacy classes, and hence 61 irreps. 56 of these have 
been identified as the irreps r " 40 " 123 of the T > mesonic rot at ion- shift group. 
The remaining 5 are obtained by inducing with the remaining non-trivial ir- 
rep of CL(4). This gives the five real irreps 80, 81, 82, 83, 16. The reductions 
of the T = baryon irreps at T > are — 

8^80, 8' ^8 2 , 16-^ 8 + 8 2 . (21) 

Under D%, we find that 

8 -> A+ + A^ + A+ + Az +E+ + E-. (22) 
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The components of the quark field which carry different representations of 
TS and TS_ are shown in Table |. 

Given the classification of baryon operators in || it is a simple matter 
to construct the D\ irreps from them. We do not present a detailed table, 
because the state of the art in measurements with baryons has not progressed 
far beyond the measurement of the purely local operators even at T = 0. 



4 Free Field Theory and Beyond 

In a free field theory of staggered Fermions, the symmetries of the Hamil- 
tonian become symmetries of the field configuration. As a result, there are 
many more degeneracies among the screening masses than in the general 
problem. The only relevant group turns out to be TS , in the sense that 
all correlators in the same irrep of TS have degenerate masses, even if they 
descend through different irreps of RF and GRF . However, the degenera- 
cies are even higher than would be predicted by the application of the group 
theory of TS . 

The local "mesons" in the Af irrep of D\ have been analysed extensively 
in free field theory (FFT) 0, [J. It is known that the correlation functions 
show the typical even-odd structure of staggered Fermion correlators. The 
screening mass, /i, in this channel is 



fia = 2Jm 2 a 2 + sm 2 [-£-) -> = 2nTa. (23) 



Here a is the lattice spacing, iV T is the lattice size in the Euclidean time 
direction, and m is the quark mass. The limit is taken for small m and large 
iV T , and is equal to twice the minimum Matsubara frequency, 2nT H (|, [OJ. 
Finite size effects are clearly strong, even in FFT, and have been analysed 
before § §. 

The remaining local "mesons" are in the Bf irrep of D%. In FFT these 
correlation functions vanish. This is easy to understand. In FFT the x and y 
direction propagators are exactly equivalent, and hence their difference (see 
Table §) cancels. 

The non-local mesons also divide into two groups. The Af , Af and E ± 
irreps of D\ give rise to screening masses according to eq. (|2~3"D . In contrast, 
the Bf and Bf irreps vanish in FFT. Deviations from any of these free field 
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theory results can be used as a measure of the interaction strength between 
Fermions. 

It is interesting to recall past measurements || f|, [J. Screening masses 
have been measured with four correlators — 

• the D\ irrep Af descending from the 1^ of the GRF , and 

• a linear combination of the two pairs of D\ irreps Af descending from 
the \\± and 2| ± of the GRF (see Table g). 

Since they all belong to the same irrep of TS_, we might expect them to give 
the same screening mass in FFT. Extensive numerical work was performed in 
the 4- flavour M and quenched || 577(3) theories at T « T c , T = 3T c /2 and 
T = 2T C . It was found that the masses within each group become degenerate 
already quite close to T c , but the masses of the two groups differed by about 
10% even at 2T C . In the light of the preceding calculations, the natural 
explanation for this observation is that there are residual interactions. 

A minor controversy persists in the interpretation of these observations. 
The two viewpoints can be summarised as — 

1. Fermions in QCD at T > T c are weakly coupled, since the Af screening 
masses coming from the and the 2^ agree very well with eq. (]23f) . 
This picture becomes a better approximation with increasing T. 

2. QCD at T c < T < 2T C is not very different from that at T < T c , since 
the Af from the 1q" ± is quite different from eq. (|2"3"D. Furthermore, the 
Af from the ± and the 2^ are just the same as the vector /pseudo- 
vector at zero temperature. 

To resolve this controversy one needs first to measure the Bf coming from 
the 2^ and check whether it agrees with either of the models above. Beyond 
this, one needs to measure screening masses in other irreps of the GRF which 
are in the same irreps of the continuum group as the measured local Fermion 
bilinears. This would check whether the second model is viable or not. 

5 Summary 

In this paper we have studied the symmetries of the transfer matrix for 
staggered Fermions. By using the method of induced representations we have 
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reproduced old, known results for the spin-flavour symmetry group and its 
irreps at zero temperature || |10[ . This method allows a simple generalisation 
to the analysis of the spatial direction transfer matrix at finite temperature. 
Our main result is the identification of the irreps of the finite temperature 
symmetry group given in Table || 

Using this, we have decomposed the T = irreps into the T > irreps. 
This reduction is shown for "mesons" in Tables ^| and [|. For quarks the 
reduction is given in Table [|. The reduction for local "baryons" may be read 
off from the same table. 

The phenomenon of "parity doubling" at high temperatures is allowed 
by the group theory. In a free fermion theory, in fact, the degeneracies are 
much higher — the screening masses are classified by the point group of the 
spatial slice of the lattice. In an interacting theory this is not true; the 
descent of each irrep through the chain of enveloping groups (see eq. [TT|) is 
important. This has been seen in the Af irrep of the point group. Whether 
the physics at T > T c is a small perturbation around the free theory or the 
zero temperature theory (or something else altogether) can be explored by 
studying other irreps. The is a good candidate because it is also built 
from local fermion bilinear operators. 

I would like to thank F. Karsch and E. Laermann for discussions when 
this work was started, and T. Venkataramanna for patient tutorials on the 
use of the theory of induced representations. 
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A Induced Representations 



Orbit 


Isotropy 


Dimension 


Multiplicity 


(0,0) 


s 3 


1 


2 






2 


1 


(0,1), (0,2), (0,3) 




3 


2 


(1,0), (2,0), (3,0) 


^2(^12) 


3 


2 


(1,1), (2, 2), (3, 3) 


ZiiRvi) 


3 


2 


(1,2), (1,3), (2,1), 








(2, 3), (3,1), (3, 2) 


{E} 


6 


1 



Table 6: Induced representations of G — (v ® V) cxi S3. A representative 
member of each orbit is underlined. The trivial isotropy group is denoted by 
{E}. This construction reproduces the result of [fTO ] . 



The method of induced representations for a semi-direct product group 



G = N x H , for Abelian N, can be found in fl4| . Here we quote the results 
required in this paper. The dual, N (set of equivalence classes of irreps of 
N), is isomorphic to N, and the action of H on N is isomorphic to its action 
on N. 

Under the action of H, the dual breaks up into disjoint orbits O iy i.e., 
N = ©jOj. Examine the isotropy group, Hi C H of one representative 
Xi € Oi. We need to know two cases — 

• When the orbit Oi has only one element, i.e., Hi = H, then the induced 
representations are precisely the irreps of H. 

• For other \i when Hi is Abelian, the dimension of the induced represen- 
tation is the number of elements in the coset H/ Hi, and the multiplicity 
of such irreps is given by the number of classes in Hi. 

All 4 irreps of the Viergruppe, V, are one-dimensional. We label them by 
the numbers 0, 1, 2 and 3. The trivial irrep is called 0. The irrep labelled 
k (7^ 0) has x{E) = xi-^k) = 1 an d the other two characters —1. Irreps of 
the direct product V <8> V are labelled by the ordered pair (k, I) where k is 
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Orbit 


Isotropy 


Dimension 


Multiplicity 


(0,0) 




1 


2 


(0,3) 


z 2 


1 


2 


(3,0) 


z 2 


1 


2 


(3,3) 


z 2 


1 


2 


(0,1), (0,2) 


{E} 


2 


1 


(1,0), (2,0) 


{E} 


2 


1 


(1,1), (2, 2) 


{E} 


2 


1 


(1,2), (2,1) 


{E} 


2 


1 


(1,3), (2, 3) 


{E} 


2 


1 


(3,1), (3, 2) 


{E} 


2 


1 



Table 7: Induced representations of G = (V <S> V) txi Z 2 . A representative 
member of each orbit is underlined. The trivial isotropy group is denoted by 
{E}. 



an irreps of the first factor and I of the second. S3 has two one-dimensional 
irreps (the trivial and the sign) and one two-dimensional irrep. Z 2 has two 
one-dimensional irreps, the trivial and the sign. These are the only inputs 
into the construction of the irreps we require. The construction of the irreps 
of G and G (see eqs. [15] and [H]), follow from the rules above, and are given 
in Tables [5] and [7] respectively. 



B Character Tables 



In this appendix we define the irreps of G and G (see eqs. [15], [H]) by writing 
down the character tables. 

In Appendix |A| we showed that the group G of eq. ([15]) has 10 irreps. 
Hence the 96 elements of the group fall into 10 conjugacy classes. For every 
g G G, we can write uniquely g = uvs, where u G V(Xk), v G V(Rl l ), 
and s G S3. We use the notation Rki and C respectively for the operators 
in 5*3 which permute the pair kl and make a cyclic shift. In terms of this 
decomposition, the 10 conjugacy classes are — the identity E, Ui (with i = 
1,2,3), Vi, u ■ v (meaning by w»Ut), UiVj (with j ^ i, and denoted uv), C, 
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G 





E 


u 


V 


u ■ V 


uv 


C 


R 


uR 


vR 


uvR 






1 


3 


3 


3 


6 


32 


12 


12 


12 


12 


1 


A 1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


V 


A 2 


1 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


2 


E 


2 


2 


2 


2 


2 


-1 














3 


Ft 


3 


3 


-1 


-1 


-1 





1 


1 


-1 


-1 


3' 


F 2 


3 


3 


-1 


-1 


-1 





-1 


-1 


1 


1 


3" 


Fx 


3 


-1 


-1 


3 


-1 





1 


-1 


-1 


1 


3"' 


F! 


3 


-1 


-1 


3 


-1 





-1 


1 


1 


-1 


3"" 


E + A 1 


3 


-1 


3 


-1 


-1 





1 


-1 


1 


-1 


3"" 


E + A 2 


3 


-1 


3 


-1 


-1 





-1 


1 


-1 


1 


6 


F 1 + F 2 


6 


-2 


-2 


-2 


2 


















Table 8: The character table for G = (V(u) <8> V(v )) M S3. The second line 
of the table gives the number of operators in each class. 

Rki, UiRij (denoted uR), ViRij (denoted vR), UiVjRij (denoted uvR). The 
character table is constructed by standard methods, and given in Table [8[ 
The reduction of irreps of G to that of the subgroup of cubic rotations O is 
performed by inspecting the characters^] of the conjugacy classes of O. These 
are the rotations by 7r/2 (R), rotations by 7r (v), rotations by 27r/3 (C) and 
the remaining 2-fold rotations (vR). 

The construction, in Appendix ^|, of the irreps of the group G of eq. 



( |17D shows that there are 14 conjugacy classes. For every g G G we have 
the unique decomposition g = uvs, where u G V(Xk), v G V(R^. l ), and 
s G Z 2 . Using the notation R for the nontrivial element of Z 2 , we can write 
the conjugacy classes as E, u = 113, v = V3, Ui (with i — 1, 2 and denoted u), 
v (in the same notation), uv, UiVi (with i — 1, 2 and denoted u • v), uv, vu, 
uv (meaning uflj with i ^ j — 1,2), R (which is equivalent to uvR, uR and 
vR), uR (also equivalent to vuR), vR (also equivalent to uvR), and uvi?. 
The character table is given in Table ^. The decomposition to irreps of -D4 
needs the identification R\ 2 = v i the 7r rotations about the x and y axes are 
v, R\ 2 = R and the remaining 2-fold symmetries are vi?. 



3 Character tables for Oh and D\ may be found in [fLJ 
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G 


£4 


E 


■u 


f 


uv 


u 


V 


uv 


f u 


u ■ V 


uv 


R 


uR 


vR 


uvi? 


1 


1 


1 


1 


2 


2 


2 


2 


2 


2 


4 


4 


4 


4 


lo 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


li 


Bi 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


1 2 




1 


1 


1 


1 


-1 


-1 


-1 


-1 


1 


1 


1 


-1 


-1 


1 


1 3 




1 


1 


1 


1 


-1 


-1 


-1 


-1 


1 


1 


-1 


1 


1 


-1 


1 4 


Ai 


1 


1 


1 


1 


-1 


1 


1 


-1 


-1 


-1 


1 


-1 


1 


-1 


1 5 




1 


1 


1 


1 


-1 


1 


1 


-1 


-1 


-1 


-1 


1 


-1 


1 


1 6 


A 2 


1 


1 


1 


1 


1 


-1 


-1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 7 


#2 


1 


1 


1 


1 


1 


-1 


-1 


1 


-1 


-1 


-1 


-1 


1 


1 


2 




2 


-2 


-2 


2 














2 


-2 














2i 


E 


2 


-2 


-2 


2 














-2 


2 














2 2 


E 


2 


2 


-2 


-2 


2 








-2 




















2 3 


E 


2 


2 


-2 


-2 


-2 








2 




















2 4 


A 1 + B 1 


2 


-2 


2 


-2 





2 


-2 























2 5 


A 2 + B 2 


2 


-2 


2 


-2 





-2 


2 
























Table 9: The character table for G = (V(u) ® V(v)) dx Z 2 . The second line 
of the table gives the number of operators in each class. 



18 



References 

[1] F. Brown et al, Phys. Rev. Lett., 65 (1990) 2491; 
Y. Iwasaki et al, Phys. Rev., D 54 (1996) 7010; 
A. Peikert et al, e-print hep-lat/901002. 

[2] B. Grossmann et al., Nucl. Phys., B 417 (1994) 289; 

S. Datta and Sourendu Gupta, Nucl. Phys., B 534 (1998) 392. 

[3] C. DeTar and J. Kogut, Phys. Rev. Lett, 59 (1987) 399; Phys. Rev., D 
36 (1987) 2828; 

S. Gottlieb, et al., Phys. Rev. Lett, 59 (1987) 1881; 
A. Gocksch, P. Rossi and U. M. Heller, Phys. Lett., B 205 (1988) 334; 
T. Hashimoto, T. Nakamura and I. O. Stamatescu, Nucl. Phys., B 400 
(1993) 267; 

Ph. de Forcrand et al., e-print |hep-lat / 990 1 1 7| . 

[4] K. Born, et al, (The MT C collaboration), Phys. Rev. Lett., 67 (1991) 
302; 

G. Boyd, S. Gupta and F. Karsch, Nucl. Phys., B 385 (1992) 481. 

[5] S. Gupta, Phys. Lett., B 288 (1992) 171. 

[6] C. Bernard, et al, Phys. Rev. Lett., 68 (1992) 2125. 

[7] T. H. Hansson and I. Zahed, Nucl. Phys., B 374 (1992) 277; 
V. Koch et al, D 46 (1992) 3169. 

[8] H. Kluberg-Stern et al, Nucl. Phys., B 220 (1983) 447. 

[9] M. F. L. Golterman and J. Smit, Nucl. Phys., B 255 (1985) 328. 

[10] M. F. L. Golterman, Nucl. Phys., B 273 (1986) 663. 

[11] R. Altmeyer et al, (The MT C collaboration), Nucl. Phys., B 389 (1993) 
445. 

[12] See for example, P. Arnold and L. G. Yaffe, Phys. Rev., D 52 (1995) 
7208. 



19 



[13] M. Hammermesh, "Group Theory and its Applications to Physical Prob- 
lems", Addison- Wesley, Reading, Massachusetts, 1962. 

[14] B. Simon, "Representations of Finite and Compact Groups", American 
Mathematical Society, 1996. 

[15] V. L. Eletskii and B. L. Ioffe, Sov. J. Nucl. Phys., 48 (1988) 384. 



20 



